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Abstract 

CAKE (Cosmic Abundances below Knee Energies) was a prototype balloon experiment for the determination of the charge spectra 
and of abundances of the primary cosmic rays with nuclear charge Z>10. It was a passive instrument made of layers of CR39® 
and Lexan®/Makrofol® nuclear track detectors; it had a geometric acceptance of ~0.7 m 2 sr for Fe nuclei. Here, the scanning and 
analysis strategies, the algorithms used for the off-line filtering and for the tracking in automated mode of the primary cosmic rays 
are presented, together with the resulting CR charge distribution and their abundances. 
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Qh The determination of the elemental abundances of cosmic 
q rays (CRs) observed in the solar system provides valuable infor- 
$-h mation on the nature of their sources, acceleration mechanisms 
c/3 and propagation in the Interstellar Medium (ISM). While there 
, ^ , is a general consensus for supernova shocks being the main 
acceleration sites of CR particles with energy per nucleon up 
CN to 10 15 eV, still unresolved is the question of which are their 
sources. 

■ From measurements of the abundances made near the Earth it 
ly^ is possible, by modelling the propagation processes in the ISM, 
fT) to infer the elemental abundances at CR sources as well as to 
place constraints on the nature of their original environments. 
^— H The present data cannot select one model between coronal ma- 
0> terial of later-type stars (which are rich of ions in a gas-phase) 
and Interstellar Medium in the proximity of Supernovae Rem- 
^ nants. Two proposed CR selection mechanisms based on the 
elemental First Ionization Potential (FIP) and on the volatil- 
rS ity/refractory state, i.e. on the condensation temperature T c , 
could justify the first scenario with respect to the second. The 
measured CR abundances show a strong prevalence of the re- 
fractory elements (which are generally characterized by low 
FIP and large T c , e.g. Mg (Z = 12), Ca (Z = 20), Fe (Z = 26) 
[1]) with respect to the volatile ones (characterized by high FIP 
and small T c , e.g. H (Z = 1), He (Z = 2), O (Z = 8)). The ap- 
parent anti-correlation between FIP and T c prevents to choose 
between the two mechanisms. Such degeneracy could be solved 
by measuring the abundances of the few volatile low-FIP ele- 
ments and the rare refractory with high-FIP ones [2]; e.g. Ga 
(Z = 31), Ge (Z = 32), Rb (Z = 37), Pb (Z = 82), Bi (Z = 83) 
[1]. 

A number of signatures exists in the abundance spectrum of 
heavy nuclei (30 < Z < 74) [1]; in order to reach confident con- 
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elusions it is necessary to achieve high-resolution charge mea- 
surements and high statistics. Elements with Z > 28 are rare 
and large exposure factors (collecting area x exposure time x 
angular acceptance) are required. For example the abundances 
of Z > 30 nuclei is 1000-=- 10000 times smaller than Fe nuclei, 
whose flux is ~1 cm" 2 sr _1 day _1 . 

Since many years high altitude balloon flights provided an al- 
ternative to more expensive (even in terms of realization time) 
satellite based experiments. Flight durations have been grow- 
ing from few days to several weeks on recent Long Duration 
Balloons (LDB) flown from Antartic regions e.g. [3], and plans 
are underway to provide missions up to months (Ultra LDB). 

Plastic passive detectors are promising since they can pro- 
vide good charge resolution, and are certainly suitable for bal- 
loon flights, as they can be used in large collecting areas and 
limited weight. The main limitation is that the detector has to 
be recovered after the flight for the analysis with microscopes. 
Such type of technique was already used succesfully in several 
experiments on balloons [4-6] and in space [7-10] 

Nuclear Track Detectors (NTDs) like CR39® can detect 
charged particles with Z/J3 > 5 with very good charge reso- 
lution [11], whereas Lexan®/Makrofol® polycarbonates have 
higher thresholds, Z//3 > 50 [12]. Due to the large flux of low 
charged particles the analysis of CR39® detectors and the de- 
termination of the relative abundances can be time consuming. 
Thus the necessity to develop an automatic system for reading 
and measuring the track characteristics. 

The purpose of the CAKE balloon experiment from the ASI 
(Agenzia Spaziale Italiana) was to test the implementation of 
new scanning and analysis techniques for future LDB flights of 
large area detectors (tens of m 2 sr) [13]. We report here about 
the scanning procedures and algorithms developed to analyse 
the CR39® sheets exposed during a trans-mediterranean flight. 
The obtained primary charge spectrum at the top of the atmo- 
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sphere is presented. 

2. The detector 

CAKE was composed of 80 multilayer stacks of CR39® 
and Lexan® sheets. Each "standard" stack was of dimension 
11.5x11.5 cm 2 , with 10 sheets of Lexan®, 0.25 mm thick, 
sandwiched between 20 layers of CR39®, 0.7 or 1.4 mm thick. 
Ten stacks were also assembled with Al or Pb foils to study 
the fragmentation of the impinging particles. Each Al or Pb 
foil was 0.5 mm thick. The thickness of a standard stack was 
~2.3 g/cm 2 while it was ~ 3.0 and 4.8 g/cm 2 for a stack with 
Al and Pb targets, respectively. Five stacks were inserted into 
an aluminium cylinder, as shown in Fig. 1, and four cylinders 
were loadged into aluminium "boxes" with an internal insula- 
tion foam coverage. The entire experimental set-up composed 
of four "boxes" was lodged on the payload of a stratospheric 
balloon. The cylinders and the boxes were designed to preserve 
the sheets in an air environment at normal preassure during the 
flight. The full exposed area was ~1 m 2 . 




Figure 1 : Composition of a cylindric aluminium container, and one as- 
sembled tray with 5 stacks. 



3. The test flight 

The balloon was launched from the Trapani-Milo base 
(12.5°E 32.92°N) of ASI and landed in central Spain after 22 
hours of flight. The plafond altitude was 37-^40 km (3-^3.5 
g cm" 2 ) and was kept for about 17.8 hours. Along the bal- 
loon trajectory the mean vertical rigidity cut-off was ~7 GV 
(limiting the minimum energy of the impinging CRs to E>2.6 
GeV/nucleon) [14]. The gondola was not azimuthally con- 
trolled and thus no East- West effect could be searched for. 

During the balloon flight the cylinders had an internal pre- 
assure of 1 Atm, and the inside temperature never exceeded 
34° C. Studies performed in a controlled heating room, deter- 
mined that CR39® keeps its sensitivity unchanged from -50°C 
up to 50°C[15]. 

4. Data Analysis 

After the flight CR39® sheets of standard stacks were etched 
in a 6N NaOH solution at (70+0. 1)°C for 30 hours; the total 



number of detector foils that were chemically etched was ~280. 
The detectors were measured by an automated SAMAICA 
(Scanning and Measuring with Automatic Image Contour Anal- 
ysis) image analyzer, i.e. an integrated DAQ system composed 
of an optical microscope controlled by a personal computer, a 
built in frame grabber, and a positioning stage [16]; the image 
analyzer was a product of ELBEK-Bildanaly se GmbH; see Fig. 
2. 




Figure 2: The used automatic measuring system. 



Chemical etching of NTDs results in the formation of conical 
pits (etch-pits) on both sides of a detector foil. The surface base 
of the etch-pits has an elliptical shape. 

For each exposed stack four CR39® detector plates were 
measured in automatic mode. For each etch-pit ("track") on 
the upper surface of a CR39® plate the system provides the 
measurement of the major and minor axis, polar angle, area of 
the base of each cone, plus its position on the foil. The scanned 
area of each foil was reduced to 8x8 cm 2 to avoid edge effects; 
a single detector sheet was measured in ~12 hours. 

By determining the CR39® bulk etching rate v# [17] it was 
possible to compute the reduced etching rate p = vt/vb, where 
vt is the etching velocity along the particle trajectory. Expo- 
sures of CR39® stacks to relativistic heavy ion beams at ac- 
celerators [18] allowed a correlation of the measured p with 
the charge number Z of the impinging ions [19]. We had at 
our disposal three CR39® calibrations for three different v# 
values, namely 1.1, 1.15, and 1.21 jim/h [20]. The interval 
1.05 jim/h < vb < 1.25 jim/h was used for accepting the bulk 
velocity of the etched foils. 236 out of 280 detectors were se- 
lected and analyzed. The parameters for the automatic mea- 
surement of the detectors were set to obtain a minimum bias in 
the track measurement; the residual noise was treated off line 
with the implementation of an appropiated filter. The scan effi- 
ciency was > 90% for the detection of tracks from nuclei with 
charges Z/J3 > 9 [15]. A filter was implemented in order to 
classify the measured tracks using a multilayer perceptron neu- 
ral network [21] trained on a data sample of -3000 cosmic ray 
tracks measured manually. The filter selection is based on two 
features of the measured tracks: the track ellipticity (defined as 
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the ratio of the minor to major semi axes of the base etch-pit) 
and its central brightness (a quantity measured in arbitrary units 
by the optical system). The filtered data set are ~ 60% of the 
raw data set [22] . 

Semi-automatic scans were made by operators to check the 
goodness of the automatic measurements and the filter; for each 
track it was verified if it was passing through both sides of the 
detector foil. The automatic scan is —10 times faster than the 
semi-automatic scan. 




Figure 3: A stack composed of 4 CR39® sheets. Charged primary CRs 
pass through the stack with linear trajectories, like the trajectory at 
the left. The tracks in consecutive detectors are matched and linked 
recursively. Selected tracks are those passing through > 3 detectors. 
Tracks with non linear trajectories, as the trajectory on the right, are 
rejected 

The trajectory of each candidate ion was followed through 3 
CR39® sheets; Fig. 3 shows different tracks crossing several 
foils of the same stack. Using tracking algorithms only tracks 
crossing the plates in linear trajectory were selected (cosmic 
ray particles with energies E > 3 GeV/nucleon are highly pen- 
etrating and scattering is negligible). A cosmic ray "event" was 
defined as a track passing through at least three consecutive de- 
tector plates as shown by the trajectory on the left in Fig. 3. 
Each candidate event was tracked individually by means of the 
fiducial area recursive method [23]. Starting with a "track" on 
the uppermost foil of the stack the method defines a limited area 
on the following foils where corresponding tracks are searched 
for. This procedure is applied recursively up to 3 etched foils on 
the stack. The fiducial area sketched in Fig. 4 was determined 
considering the experimental uncertainties of the track inclina- 
tion (in Fig. 4 6 is the complementary of the zenith angle) and 
its orientation (cp is the azimuth angle) with respect to the foil 
uppermost surface. All the uncertainties of the measured vari- 
ables are considered, i.e. the etch-pit base area (±10 yum 2 ), the 
foil thickness (at most ±50 jim) and the mis-alignment among 
the foils in the stack. The size of the fiducial area was <1 cm 2 . 
Details about the tracking and selection procedures are given in 
ref. [15,22]. 

5. Results 

In Fig. 5 is shown the unfolded charge distribution of the 
CR events detected in CR39®. The range of the reconstructed 
charge number Zis5<Z<30. In the same figure are shown 
the Gaussian distributions of each element (represented with 
thin black curves) which all added together fits the experimental 





"t"\ p eosAO i e \~2 

If 


'. : ™ 

tanG 

z AO 

/'^£\ sin 2 


V 

(7) 



Figure 4: "Fiducial" area where a candidate CR track is searched for in 
subsequent foils. The area takes into account the experimental uncer- 
tainties in zenith (90° -0), azimuth <p, the distance between consecutive 
foils z, and the mis-alignment beween them e. 



histogram (represented with the thick black curve). The Gaus- 
sian parameters (height and width) were computed by means 
of a global fit of the sum of 25 Gaussians to the experimental 
distribution. The fit was computed using the least squares min- 
imization package [24]. The standard deviations of the fitted 
curves were in the range 0.4 0.6 Z. By integrating over each 
individual Gaussian distribution we obtained the distribution of 
observed events at the detector level shown in Fig. 6. 
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Figure 5: Experimental charge distribution (histogram). A global fit 
(thick black curve) is applied using a Gaussian approximation for the 
individual peaks (represented by the thin curves). 



5.7. Fragmentation processes 

In order to compute the CR charge distribution at the top of 
the atmosphere, fragmentation processes have to be considered 
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[25]. Fragmentation of cosmic ray nuclei may occur in the in- 
teractions with the residual atmosphere; the mean atmospheric 
thickness at the CAKE trajectory was 3.25 glcm 2 . Fragmenta- 
tion also occurs in the interaction with the aluminium contain- 
ers housing the detector stacks, shown in Fig. 1; the mean Al 
thickness was 2.7 glcm 1 . Fragmentation on air and on Al were 
modelled with Monte Carlo simulations using Geant4 [26] . The 
model takes into account the interactions of CR nuclei with the 
residual atmosphere and the detector shields. Primary CR en- 
ergies were sampled assuming the energy dependence of the 
flux O = E~ 2J [27] in the range 5 < Z < 28. With the sim- 
ulations the fragmentation coefficients of heavy elements into 
lighter species were obtained, calculated as in ref. [28] and 
[29] ; then with the fragmentation matrix the charge distribution 
at the top of the atmosphere was obtained. In Fig. 6 the black 
histogram shows the charge distribution at the detector level; 
the solid red histogram is the charge distribution at the top of 
the atmosphere. 
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Figure 7: Angular acceptance multiplied by the exposure area and the 
exponential attenuation factors as a function of the charge Z. 
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Figure 6: Charge distribution of the unfolded distribution with a resolu- 
tion of \e in the range 5 < Z < 30. The dotted black histogram is the 
distribution of events for each species at the detector level. The solid 
red histogram is the charge distribution at the top of the atmosphere 
obtained by correcting for fragmentation processes in the aluminium 
containers of the detectors and in the residual atmosphere. 

5.2. Angular Acceptance 

The geometrical acceptance of the C39 detectors varies with 
the charge Z [19] as there is a limiting angle 6c for the formation 
and observation of a track in the detectors 



sin(0 c) = -= V J- 
P vt 



(1) 



By including the exponential attenuation of CRs (i.e. by consid- 
ering the residual atmosphere, the Al containers of the detectors 
and the insulating foam which filled the vessels) the geometri- 
cal factor of the detector varies as shown in Fig. 7. 

5.3. Relative Abundances at the top of the Atmosphere 

The abundances of nuclei at the top of the atmosphere, were 
evaluated considering the angular acceptance and correcting for 



the flux attenuation in air and in the detector boxes. The re- 
sulting values are shown in Fig. 8 together with the measure- 
ment reported in ref. [30] at a similar energy threshold (2.6 
GeV/nucleon). The agreement is satisfactory for 12 < Z < 26. 
For Z < 12 there is a disagreement due to the filter and algo- 
rithm rejecting small tracks affecting mainly low charged parti- 
cles. The differences at the left of the Fe peak are probably due 
to a partially unsatisfactory correction for fragmentation. At the 
right of the Fe peak the number of events is too small to even 
consider a significant result. 

In Fig. 8 error bars in CAKE's data represent systematic plus 
statistical standard deviations added in quadrature; the error 
around the Iron relative abundance is ~ 20%. 
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Figure 8: Measured CAKE abundances relative to iron (set equal to 
100), as a function of the charge (square points and red solid line). 
Even charge elements are quoted in the plot. Error bars correspond 
to systematic and statistical standard deviations added in quadrature. 
Experimental measurements of a similar experiment compiled in ref. 
[30] are also plotted as black points interpoled by a dashed black line. 
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6. Conclusions 

This paper was focused on the development of automatic 
analysis procedures for determining the chemical composition 
of cosmic ray nuclei using nuclear track detectors. The NTD 
technique offers good charge resolution and is attractive for 
the exploration/determination of the charge spectrum of cosmic 
particles at reasonable cost. The explotation of this technique 
in large area experiments relies on automatic image analyzers 
for scanning and measuring the detector foils as well as on the 
availability of LDB and ULDB flights. 

By making use of a short stratospheric balloon flight we had 
the opportunity to study and develop several tools to obtain a 
charge spectrum of cosmic ray particles with Zip > 10 at the 
top of the atmosphere. Our measurement is in fair agreement 
with other observations. 
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